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We  shall  investigate  the  stability  of  a  solution  of  the  following 
system  of  differential  equations  which  arose  from  the  investigation  of  shocked 
waves  by  Professor  Taub: 
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Here  and  in  the  sequel  we  shall  use  the  same  notations  as  in  a  paper  "Deter- 
mination of  Flows  Behind  Stationary  and  Pseudo-stationary  Shocks"  by  Professor 
Taub. 

To  discuss  the  problem  the  most  important  coefficients  in  this  system 
are  A.  .   (i,j  =  1,2)  and  they  are  given  by  (5„ll): 
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The   characteristic  roots   of  this  matrix  are   the   solutions   of 


+2  M         2   sin  u  1 .  n 

Wn  2        2    x      x  *     .2-  2        2    »  U' 

\(1  -   m  cos   mJ  a.    (1  -  m  cos   uj 


namely 


sin  u cos   u  ■/  2  ,    . 

a+  ■-  ~Z 2~T7  i   — 2— TT  V»   -  1  (3) 

—         \(1  -  m  cos  |i)  \(1  -  m  cos  p.) 


-1- 


We  may  write  our  system  of  equations  (l)  in  the  following  form: 
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where  a. 4  and  b.  are  known  functions  of  x, ,  ....  xc 
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Our  case  will  be 
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The  characteristics  values  of  the  matrix  (a.  .)  will  be  Ct     and  four  0's. 
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Suppose  that  there  are  variations  Ox.  of  the  solutions.   Then  these 
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We  assume  <$x.  takes  the  form 
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where  a   is  real.   Our  system  of  equations  (h)   will  be  stable  if  the  imaginary 
part  of  £   is  <  0.  We  shall  consider  each  x.  to  be  fixed,  so  that  each  co- 
efficient a.  .  (and  b.)  is  a  constant,  and  a  is  sufficiently  large,  so  that 
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Substituting  (7)  into  (6),  we  get  a  system  of  linear  equations  between  a  andx7 

with  constant  coefficients:   namely 

We  have  used  an  approximation  that  dx,  are  so  small  that  their  square  may  be 
neglected.  From  our  convention  a   is  so  large,  and  hence  (8)  may  become 
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Hence  CK^  must  be  one  of  characteristic  values  of  the  matrix  (a. .),  i„e„ 

+ 
In  our  case 


A=   <i\(sin  (J.  4  cos  M-Vm  -  1  )  .  (9') 

Hence  the  system  (l)  _is  stable  if  m  ^  1. 
In  case  the  flow  is  subsonic  the  investigation  failed  to  produce  any  definite 
statement  on  the  stability  of  equations. 

In  the  following  we  shall  discuss  the  similar  problem  for  the 
difference  equations  derived  from  (l).  Again  we  shall  take  our  equations 
in  the  form  (h) . 
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Replace  = —  and  -5 —  by  the  difference  quotients 
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and  use  the  notations  such  as 

xi'^  =  xi(n4u,i^v). 
Our  basic  difference  equations  will  be 
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a  is  real  and  is  to  be  assumed  sufficiently  large.   The  equations  (10)  will  be 
stable  if  |  J  I  ^  1  for  all  (large)  values  of  a.   Substituting  (il)  into  (10)  we 
get 
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The  principal  term  of  the  right  hand  side  is  the  first  one  as^du,^Jv->  0  keeping 
<du/dv  in  a  finite  value.   Hence  again  using  (ll) 
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This  relation  may  be  obtained  from  (8')  since  (8')  gives  an  approximation  of 
variations  when  a   is  large.   The  relation  (12)  shows  that  —r~  is  a 


characteristic  value  of  the  matrix  (a.  .) 
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ition  for  the  stability  is  £  <.  1  for  all  s  with  s  =  1 


Now 


5  =  1  h —  t^(s  -  l),  so  that  the  above  condition  is  satisfied  if  and  only  if 
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The  positiveness  of  the  right  hand  side  is  guaranteed  for  a  shock  of  non-zero 
strength,  since  in  that  case  we  have 
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If  the  flow  is  subsonic:   m  <  1,  then  the  system  (10)  of  difference  equations  is 
unstable,  i.e.  the  solution  does  not  converge  to  a  solution  of  our  original 
equations  (l)  whenZlu,  Av  — *  0. 

There  are  many  ways  to  set  up  the  difference  equations  corresponding 
differential  equations.   Integrating  both  sides  of  (k)   with  respect  to  u  and 
applying  the  trapezoid  rule  to  the  right  hand  side,  we  get 
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where  x.   =  x.(n£u,v),     b.    =  b.(x, ,    ...,   x,  )   etc.      Considering  cfx.    of  the   form: 
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dx.  =  ox.s  e    where  s  =  e  '    ,  we  get  a  difference  equation  for  <£x.  such  as 
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In  deriving  (15)  we  assumed  that  dx.  are  so  small  that  their  square  may  be  neglected, 
The  stability  condition  for  the  system  (l^+)  is  that  the  real  part  of^  is  less  than 
or  equal  to  0  whenever  |  s|  =1.   In  the  actual  form  of  (15)  it  seems  difficult  to 
give  an  explicit  condition,  so  we  assume  here  that^u  has  been  taken  so  small  that 
we  may  replace  (15)  by 
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If  this  approximation  is  permissible,  then  we  can  solve  A   as  a  function  of  s: 
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where  a  are  characteristic  values  of  the  matrix  (a.  . ).   Since 
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Thus  the  real  part  of^  ^  0  for  all  |  s|   =  1,  if  and  only  if  a     are  real,  i.e. 
m  >.l. 
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Summary 

The  system  (l)  of  equations  is  stable  if  the  flow  is  supersonic. 
I  do  not  know  whether  stable  or  not  in  the  other  case.   The  similar  stability 
conditions  for  the  corresponding  difference  equations  are  examined  for  two 
kinds  of  systems.   In  both  cases  the  systems  are  stable  if  the  flow  is  super- 
sonic, but  unstable  in  case  of  subsonic  flows. 
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